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Abstract 


Measurements of limb scattered radiance spectra performed with the SCanning Imaging Absorption SpectroMeter for Atmo- 
spheric CHartographY (SCIAMACHY) aboard ENVISAT are used to detect noctilucent clouds (NLCs) in the northern hemi- 
sphere during summer 2002. The NLCs are easily identified by enhanced limb radiances around 80 km tangent height. Using limb 
radiances in the UV-B/C spectral range, where multiple scattering is negligible, the spectral signature of the NLC backscatter can be 
determined directly. Mie-calculations are used to estimate particle sizes of these clouds from SCIAMACHY measurements. In 
particular, we use the information contained in the spectral exponent of the NLC limb backscatter signal. The derived particle radii 
assuming a log-normal NLC particle size distribution are 40-50 nm and they are in excellent agreement with previous particle size 


estimates. 
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1. Introduction 


The summer mesosphere still receives a significant 
amount of scientific interest, particularly related to 
noctilucent clouds. NLCs occur near 82 km altitude 
during about a 12 week period centered 2 weeks after 
summer solstice in the summer hemisphere, poleward of 
about 55° (Thomas and Oliviero, 1989), where temper- 
atures frequently drop below 150 K. The particles pre- 
dominantly consist of water ice (Hervig et al., 2001) with 
radii ranging from about 10 nm to about 80 nm (e.g., 
Witt et al., 1976; Rusch et al., 1991; Gumbel et al., 1999; 
von Cossart et al., 1999). Larger radii (up to 350 nm) 
have also been reported (Thomas, 1991). The study of 
NLCs with UV/vis optical methods has a history of 
several decades. NLCs were detected with limb viewing 


* Corresponding author. Tel.: +49-421-218-4582; fax: +49-421-218- 
4555. 
E-mail address: csavigny@iup.physik.uni-bremen.de 
(C. von Savigny). 


satellite spectrometers (OGO 4 (Donahue et al., 1972), 
UVS on SME (Thomas and Oliviero, 1989), OSIRIS on 
Odin (e.g., Petelina et al., 2001)) and imagers (WINDIT 
on UARS (e.g., Evans et al., 1995)), rocket-borne pho- 
tometers and polarimeters (e.g., Witt et al., 1960; 
Gumbel et al., 2001) and ground-based LIDARs (e.g., 
von Cossart et al., 1999; Alpers et al., 2000). The esti- 
mation of NLC particle sizes is possible from measure- 
ments of color ratios of the NLC backscatter signal 
(e.g., Rusch et al., 1991; von Cossart et al., 1999), the 
NLC extinction spectrum in solar occultation (Debres- 
tian et al., 1997) and from measurements of the NLC 
scattering phase function with photometers on spinning 
rockets (Witt et al., 1976; Gumbel et al., 2001). Gumbel 
and Witt (2001) showed in a case study that the particle 
size estimation by color ratios and measurement of the 
phase function yield consistent results. Recently, the 
determination of both radius and width of an assumed 
log-normal particle size distribution became possible 
with multi-wavelength LIDAR measurements. In this 
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context, von Cossart et al. (1999) studied 11 NLC events 
over Andoya (Norway) with a 3-wavelength LIDAR 
system. They assumed that the NLC particles are water 
ice spheres, having a log-normal particle size distribu- 
tion (PSD) 

f(r) 


In? (r/ro) 


21n? o (1) 


a V2nrlno 


where N is the number density of NLC particles and r is 
the radius of a particle. Using Mie-theory, von Cossart 
et al. (1999) found values of 7 = 51421 nm and 
æ = 1.42+0.22 for the 11 events. They also estimated 
the particle number density and obtained a value of 
N = 82 +52 cm’? based on their assumptions. To what 
extent the non-sphericity of the NLC particles (Baum- 
garten et al., 2002) and their refractive index influence 
the retrieval results remains unclear. 

Each NLC observation method has its advantages 
and disadvantages: LIDARs and rocket photometers 
generally have a much higher vertical resolution (usually 
better than about 200 m) than limb viewing satellite- 
based instruments (typically 1-3 km). Rocket-borne in- 
struments provide measurements at a single point in 
space and time only, whereas LIDAR measurements can 
in principle provide an extended temporal record of 
NLC presence, altitude, and particle size estimate for a 
given location. The advantage of satellite observations is 
their ability to provide a global picture of NLC occur- 
rence and particle size estimates. The exact altitude of 
NLCs cannot be determined with limb scattered sunlight 
observations, only an upper limit can be established. 
Moreover, observations with sun-synchronous satellites 
cannot provide information on the local solar time de- 
pendence of atmospheric phenomena, since air volumes 
at a certain latitude are always sampled at the same local 
solar time. 

The purpose of this paper is to present preliminary 
results of NLC observations and particle size estimates 
from SCIAMACHY limb measurements and to show 
that SCIAMACHY with its global coverage and its wide 
spectral range has the capability to significantly improve 
our understanding of the summer mesosphere. 


2. Brief instrumental description 


SCIAMACHY is one of 10 scientific instruments on 
board the European Space Agency’s ENVISAT satellite. 
ENVISAT was launched on March Ist 2002 from Ku- 
rou (French Guyana) into a sun-synchronous polar or- 
bit with an inclination of 97.8° and a descending node at 
10:30 local solar time. SCIAMACHY is the first in- 
strument capable of performing spectroscopic observa- 
tions of the Earth’s atmosphere in three observation 
modes: nadir viewing mode, limb viewing mode, as well 
as the solar and lunar occultation modes. For this study 


only SCIAMACHY limb observations were employed. 
SCIAMACHY has eight channels each of which con- 
tains a grating. The spectral range between about 220 
and 2380 nm is covered with a spectral resolution of 
about 0.2 nm in the UV-B/C spectral range relevant to 
this study. In the limb observation mode the Earth’s 
limb is stepped through for a tangent height range be- 
tween —3 and 100 km in steps of about 3.3 km. After 
each tangent height step an azimuthal scan is performed 
covering a distance of about 960 km at the tangent 
point. The entire atmosphere from the troposphere up 
to the mesopause is observed all the time. Thus, no 
dedicated summer mesosphere campaign is necessary for 
the observation of NLCs. SCIAMACHY’s instanta- 
neous field of view spans 2.6 km vertically and 110 km 
horizontally at the tangent point. Further information 
on the instrument, measurement modes and mission 
goals was given by Bovensmann et al. (1999). 


3. Methodology 


For the NLC studies presented here SCIAMACHY 
limb radiance observations in the UV spectral range 
(channel 1, covering the range from about 220-310 nm) 
were employed. The use of UV wavelengths is advan- 
tageous compared to visible wavelengths, and this is 
why the UV spectral range has been intensively used for 
NLC studies (e.g., Witt et al., 1976; Gumbel and Witt, 
2001; Gumbel et al., 2001). Due to the strong absorption 
within the Hartley and Huggins bands of ozone the 
detected limb scattered photons have almost exclusively 
experienced only one scattering event (e.g., Kaiser and 
Burrows, 2002). The single scattering approximation 
(Kokhanovsky, 2001) is therefore sufficient to interpret 
the measurements. This also implies that the ground 
albedo has no impact on the observed radiances. 
Therefore, the spectral signature of the NLC backscatter 
can be determined with high accuracy. 

Fig. 1 shows examples of SCIAMACHY limb radi- 
ance profiles at several UV wavelengths with NLCs 
within the field of view (FOV) (panel a), and without 
NLCs in the FOV (panel b). In the absence of NLCs the 
limb radiance increases with decreasing tangent height 
as long as one is within the optically thin regime, i.e., the 
line of sight (LOS) optical depth t is small. The observed 
limb radiance consists of different contributions in the 
presence of NLCs: (a) the molecular Rayleigh back- 
ground; (b) backscatter signal from the NLC. The limb 
radiance profiles in Fig. 1b exhibit a maximum (the so- 
called “knee”) between about 50 and 70 km depending 
on the wavelength. This is indicative of a LOS optical 
depth exceeding a certain threshold (t >> 1). The atmo- 
sphere becomes optically thick, i.e., the transmission 
from the tangent point to the observer is close to zero. 
The knee altitude is only weakly dependent on geo-lo- 
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Fig. 1. Examples of UV limb radiance profiles observed with SCIAMACHY: (a) with NLCs; (b) without NLCs in the field of view. 


cation and geophysical parameters and has therefore 
also been employed for attitude determination of limb 
sensors (e.g., Sioris et al., 2001). 

For this study Level 0 data were used that have been 
corrected for dark current and responsivity. The in- 
strument’s polarization sensitivity is not taken into ac- 
count, but this has a negligible effect on the limb 
radiance profiles, since the degree of linear polarization 
is almost independent of tangent height, as confirmed 
with the spherical vector radiative transfer model of 
McLinden et al. (2002). 

Fig. 2 shows a polar plot, of all NLCs detected with 
SCIAMACHY on August 11, 2002. Note that the data 
of only 10 of the 14 orbits performed on this particular 
day were available and are shown in Fig. 2. The rect- 
angles correspond to the area covered by the tangent 
point ground track for each limb scan. The presence of 
NLCs is indicated by solid black rectangles, hatched 


NLCs detected by SCIAMACHY on Aug. 11, 2002 


Fig. 2. NLCs detected by SCIAMACHY on August 11, 2002 are de- 
picted by solid black rectangles. Hatched rectangles correspond to 
SCIAMACHY limb scans where no NLC was detected. 


rectangles are NLC free. Obviously, NLCs only cover 
the highest latitudes on August 11, typical for the re- 
traction of NLCs to higher latitudes at the end of the 
NLC season. 

The NLCs are optically thin with visible slant peak 
optical depths of typically less than 5 x 107? (Debrestian 
et al., 1997). The total UV-Vis LOS optical depth from 
the tangent point at a tangent height of 82 km (the 
typical altitude of NLCs) to an observer in space is 
smaller than 5 x 10-7. This is valid even in the center of 
the O; Hartley bands near 250 nm. In the single scat- 
tering approximation the contribution of the observed 
limb radiance Jnic(/, @) originating from scattering by 
the NLC can be calculated in the following way (Ko- 
khanovsky, 2001) 


Inuc(A, O) = A x (A, @) x S(A), (2) 


where q(A,@) is the NLC differential scattering coef- 
ficient, © denotes the scattering angle, S(A) is the 
solar spectrum, and A corresponds to a spectrally 
neutral proportionality constant. Employing longer, 
i.e., visible wavelengths is in principle desirable, since 
then both parameters of an assumed log-normal par- 
ticle size distribution (see Eq. (1)) can in principle be 
determined (Alpers et al., 2000; von Cossart et al., 
1999). Yet, the multiple scattering source function has 
to be taken into account and this will generally not be 
known with sufficient accuracy, since surface albedo 
and cloud cover would have to be known and mod- 
eled exactly. 

The NLC sun-normalized limb scatter spectrum 
i(A, O) = ee which is proportional to g(/, O). was 
derived (see Ha, (2)) at wavelengths in the UV-B and 
UV-C (265-300 nm). 265 nm was chosen as the lower 
boundary to avoid contamination by the stronger NO y- 
bands at shorter wavelengths. The Rayleigh background 
was determined from model calculations with the radi- 
ative transfer model SCIARAYS (Kaiser, 2001). The 
simulated limb radiance profiles were scaled to match 
the observed limb radiance profiles at the “knee”. The 
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Limb radiance spectrum at 83 km 
tangent height Scan 22, Orbit 1778 .. . 
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Solar irradiance spectrum measured with 
SCIAMACHY [Skupin et al., 2002] 
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Fig. 3. Limb radiance spectrum for an NLC observation and solar irradiance spectrum (panel a) used to determine the spectral dependence of the 
sun-normalized NLC spectrum i(4,@) (panel b). Panel (b) also shows power law functions (crosses) fitted to the measured sun-normalized NLC 


spectrum (solid lines). 


Modelled spectral exponent a for a log- 
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Fig. 4. Panel (a): theoretical dependence of the spectral exponent « on the log-normal PSD parameters rọ and o within the 265-300 nm spectral 
window for scattering angle of State 23. Panel (b): color ratios c(A) obtained from measurements (solid lines) and Mie calculations (crosses) for 
monodisperse NLCs (scaled by a factor of 2 for scan 23 and 3 for scan 22) to allow a better discrimination of the different spectra). 


Skupin et al. (2003) solar irradiance spectrum (see 
Fig. 3(a)) as measured by SCIAMACHY was used as 
S(2). We found that the NLC differential scattering 
coefficient can be approximated by a power law 
q(4, 0) x 2~* with the spectral exponent «. 

Fig. 3(b) shows a log—log plot of the sun-normalized 
NLC spectrum as a function of wavelength. The abso- 
lute of the slope corresponds to the spectral exponent «. 
The spectral points nicely follow a power law. Due to 
SCIAMACHY’s hyperspectral capability the spectral 
exponents « can be determined with high fidelity. Note 
that the sun-normalized spectra i(2, ©) were smoothed 
with a 2 nm boxcar to remove the high frequency 
structures that arise mainly due to small wavelength 
misalignments. This high frequency component is not 
relevant for the determination of the spectral exponents. 
To infer information on the NLC particle sizes, the 
spectral exponents « obtained from SCIAMACHY ob- 
servations are related to particle sizes by model calcu- 
lations. For this we used the same assumptions as von 
Cossart et al. (1999) (see Section 1) and assumed that 


Eq. (1) is valid. The refractive index for ice was taken 
from Warren (1984). Fig. 4 shows model results of the 
spectral dependence of the differential scattering cross 
section for polydisperse NLCs (panel a) and monodis- 
perse NLCs (panel b). Fig. 4a is a contour plot of the 
NLC spectral exponent for the 265-300 nm wavelength 
range as a function of rọ and o for the viewing geometry 
of limb scan 23 of orbit 1778. Fig. 4b gives the color 
ratios c(A) = ay for scans 22, 23, and 25 obtained 
from SCIAMACHY measurements. Also shown are 
Mie calculations for monodisperse spheres (i.e., ¢ — 0 


in Eq. (1)). 


4. Discussion 


The spectral exponent in Fig. 4(a) apparently depends 
on both rọ and o of the assumed log-normal distribu- 
tion. A unique determination of both size distribution 
parameters requires at least three significantly different 
wavelengths (e.g., von Cossart et al., 1999). Due to 
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SCIAMACHY’s wide spectral coverage, the determi- 
nation of both ro and ø is in principle possible. Yet, this 
is difficult in practice since limb radiance profiles at 
visible wavelengths have to be employed. Generally, the 
unknown multiple scattering source function will make 
the accurate extraction of the NLC spectral signature 
difficult. Therefore, we limit ourselves to UV wave- 
lengths here. The determination of both size parameters 
is presently under investigation. If we assume o = 1.4 
(von Cossart et al., 1999) (see Section 1), then the esti- 
mated value of rọ of the NLC signature observed by 
SCIAMACHY is about rọ = 45 nm for scan 22, i.e., 
within the typical range of NLC particle sizes. For scans 
23 and 25 we obtain rọ = 50 nm and rọ = 40 nm, re- 
spectively. If we allow the width parameter o to fall 
within the range 1.2-1.6. then the inferred radii are 25- 
65 nm for scan 22, 35-75 nm for scan 23, and 25-55 nm 
for scan 25. 

If a monodisperse NLC model is used (Fig. 4(b)) 
then the radii that fit the observed normalized back- 
scatter spectra best are 92 nm for scan 22, 95 nm for 
scan 23 and 80 nm for scan 25. Interestingly, the radii 
obtained from the monodisperse model are approxi- 
mately twice as large as the radii from the polydis- 
perse model assuming a log-normal distribution and 
c= 1.4. The fact that the monodisperse radii are 
greater than the polydisperse ones is not surprising, 
since the scattering cross section scales as the squared 
volume of the particle in the Rayleigh limit, i.e., larger 
particles scatter more efficiently, and therefore they 
dominate the total scattering cross section. We found 
that the observed NLC backscatter spectra can be 
modeled very accurately with both a monodisperse 
and a polydisperse model. A similar result was re- 
cently obtained by Gumbel and Witt (2001) who 
showed that the observed NLC phase functions can 
be explained using both monodisperse and _ polydis- 
perse models. This amounts to an ambiguity of the 
method to determine NLC particle sizes and the 
question arises what model is to be used. Both models 
do not exactly represent the actual PSD. The log- 
normal PSD — which is the most common PSD used 
to infer NLC particle sizes from optical measurements 
— most likely overestimates the number of larger 
particles, since large particles actually fall out faster 
than smaller ones (e.g., Gumbel and Witt, 2001), al- 
though it nicely fits the multi-color LIDAR measure- 
ments by von Cossart et al. (1999). The monodisperse 
model is obviously inadequate. Yet, one could argue, 
— particularly since the actual PSD is unknown that 
one should employ the simplest model -that accurately 
describes the observations. Once a more realistic NLC 
PSD is found, the particle radii based on the mono- 
disperse assumption can numerically and perhaps an- 
alytically be transformed to the more appropriate size 
parameter. 


5. Conclusion 


First results of SCIAMACHY observations of 
NLCs in the northern hemisphere during summer 
2002 and estimations of the radii of the NLC particles 
are presented. The NLC backscatter spectrum i(/, 0) 
and its spectral exponent « between 265 and 300 nm 
were determined for several cases. NLC particle radii 
were estimated assuming both monodisperse and 
polydisperse ice spheres, and the observations are re- 
produced very accurately with both models. The de- 
duced NLC particle radii fall within the range 40-50 
nm assuming a log-normal particle size distribution 
with o = 1.4 and between 80 and 95 nm based on the 
monodisperse model. SCIAMACHY continues the 
long-term observations of global NLC coverage star- 
ted by SME and complemented by WINDII, and it 
has the unique potential of providing a long term 
archive of NLC occurrence, particle sizes and the 
variations thereof. 
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